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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 1333 


ON ROTATIONAL CONICAL FLOW* 
By Carlo Ferrari 

SUMMARY 


The author determines some general properties of isoenergetic rota- 
tional conical fields. For such fields, provided the physical parameters 
of the fluid flow are known on a conical reference surface E, it being 
understood that they satisfy certain imposed conditions, it is shown how 
to construct the hodographs in the various meridional semiplanes, as 
the envelope of either the tangents to the hodographs or of the osculatory 
circles . 


ANALYSIS 


1. A method for determining the field of flow about a cone of revo- 
lution, the axis of which is aligned with the direction of the impinging 
supersonic stream, which is taken to have a uniform velocity distribution 
sufficiently far ahead of the body before the conical field is created, 
was developed by Busemann (reference l) at an early date. Several years 
later the present author (reference 2) extended Busemann's procedure to 
cover the case of a cone of any shape whatsoever situated in the flow, 

so that its axis was at any arbitrary finite angle of attack; that inves- 
tigation was confined, however, solely to irrotational conditions. With 
proper alterations, nevertheless, this treatment of the yawed arbitrary- 
shaped conical surface can be applied, to the case of rotational flows . 

The purpose of the present investigation is just to give the relationships 
which permit one to draw the hodographs for the flow, in reference to 
the various meridional planes through the body axis, in the case where 
the motion is defined by a rotational field of conical flow with arbitrary 
specification of the cone shape. 

Y 

2. Upon the mere assumption of isoenergetic flow in a perfect fluid, 
one may write the equations of motion in the form (reference 3): 

^Original Italian Report appeared as Sui Moti Conici Rotazionali In. 
Onore di Modesto Panetti , published by L' Aerotecnica, Associazione Tecnica 
Automobile, and La Termotecnica, Turin, Italy, November 25, 1950- 
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-> n 2 

curl V x V = ±- 
yR 


grad S 


div 



( 1 ) 


wherein R is the universal gas constant, y is the adiabatic exponent, 
Vi is the limiting velocity obtained when the flow expands to a vacuum, 

C is the velocity of sound, S is the entropy, and V is the fluid 
velocity. 

Upon employment of a spherical coordinate system (r,0,cp), as depicted 

curl V 

in figure 1, the components of — — are taken as w r , and w q‘, 

where w r is the radial component, w f p is the component lying in the 
meridional plane and normal to the radius vector, while w q is the com- 
ponent that is perpendicular to the meridional plane. Likewise, the 
corresponding components of are denoted by v r , Vq>, and v 0 . 

Between these components there subsist the following relationships 


1 

sin cp 


dv<p 

S0 


S(v 0 sin cp) 

Sep 


= rw r 


S( rvg) 
Sr 


1 Sv r 
sin cp S© 


rwcp 


( 2 ) 


Sv r 


Sep 


S(rvcp) 

Sr 


rw 0 


Based on the assumption that the flow is conical, the following 
scalar equations are derived in a straightforward way from the first of 
the equations ( l) : 
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w eV ' w cp v 0 = 0 


- v 9 w r 


w cp v r = 


£ 1 Ss 

?R r sin cp 5 e 


( 3 ) 


v 0 w r 


v r w 0 


1. £ 5s 
r 7R dcp 


wherein £ = £/v£. 

The set of equations (3) are not independent of each other as is 
evident from consideration of equation (l) directly, but they are inter- 
related through the expression 

V x grad S = 0 ( 4 ) 

which, for a conical field, becomes 

1 dS , ) 

Sep \ sin cp 5e 

From the second of equations ( l) , upon use of the hypothesis that 
the flow is conical, and by taking into account the relationships expressed 
by equations ( 2 ) and (4'), one then obtains that 




2 v cp v 0 c* v cp 
sin cp c 2 ^0 


cot cp Vq) 





1 ^e\ 

sin cp 50 ) 


Vp &3 v cp v r 

7R dp c 2 


(5) 


This expression differs only by the presence of the rotationality terms 
from the analogous relationship derived in reference 2 previously mentioned. 
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3. By means of equations (3)> (b) , and (V) one may deduce some 
interesting properties of conical fields. Let it be assumed that one of 
the stream surfaces of the flow is conical (this will be the case for a 
field of flow arising by the action of a uniform supersonic stream 
impinging from any direction whatsoever upon a conical-shaped obstacle); 
it will be convenient to designate this surface as E c . Let the versor 

of any arbitrary general one of the generatrices of the conical sur- 
face Z c be denoted by r, then grad S x ? = 0. On the other hand, if 
the versor of the tangent to any streamline whatsoever that is traced 
upon the surface of the cone L c is' denoted by t then it is true, in 

addition, that grad S x t = 0. It is evident, therefore, that grad S 
is perpendicular to the surface Z c ; that is, the above -described conical 

surface is a surface of constant ■ entropy. 

Besides, let it be assumed that the conical flow is symmetric with 
respect to the meridional plane 0 = ± 90 ° (this will be the case already 
mentioned for the field of flow about a conical-shaped obstacle) . At 
all the points of this plane it is true that vg = 0. One then deduces, 
upon the basis of the first of equations ( 3 ) , that vg = 0 provided 
that v^ is not' zero everywhere. On account of this, and through 

utilizing the third of equations ( 3 ) it follows that — = 0. Thus even 

dp 

the meridional plane of symmetry for the conical field is itself a con- 
stant entropy surface, and at this plane the flow is irrotational as is 
easily deduced upon taking cognizance of equations (2) . 

In conjunction with the result obtained above one can derive from 
this latter fact that, for the case of flow about a conical obstacle, the 
shock wave in the two semiplanes 0 = 90 ° and- 9 - - 90 ° must produce 
the same change in direction of the stream velocities; and so the tangents 
to the trace of the shock wave in these semiplanes are symmetrically 
inclined with respect to the undisturbed stream velocity vector. Now 
let us consider an obstacle in the form of a right circular cone. Upon 

the surface of this cone it is true that — = vm = 0, and therefore one 

<30 

gets that wjp = 0. Thus the following relationship results 


1 



v 0 = 


sin cp d0 


( 6 ) 
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On the cone one can always express the v r values as a periodic 
function of 9 , and thus v r = B 0 + Z B n sin n 9 . It follows that on 
the cone's surface the peripherial velocity component is given, by 


v q - — i — Z B n n cos n0 
sin cp 

just as in the case of irrotational flow. 

Now, if we let the angle of incidence of the axis of the cone be 
denoted by (3 , then the expression for v r becomes simply 


v r = B 0 + B-^ sin 9 

if only terms of the order of magnitude of 3 are taken into account. 
Since this is true, then because B^ is proportional to 3> the B n 

p 

coefficients have to be at least as small as 3 * 

The relationship given by equation (6) may be generalized for the 
case of a cone of any shape whatsoever. It is assumed for this purpose 
that the cone's surface is divided up by a network of orthogonal coordi- 
nate lines orp and a 2 (r = const, and cp = const., respectively). 

The former of which are the intersections of the spheres with radius r 
upon the cone under consideration, while the latter are the generatrices 
of the cone itself. At an arbitrary general point P on the cone the 
length of the linear element dcr^ can ^e written as: da^ = rh^ (0) d0 

while the length of the linear element dcr 2 along the line a 2 is 
given by: dcr 2 = dr. 


The component, in the direction of the normal to the cone at the 
point P, of the curl is 


V n = & ( V l h l r ) 


dv 2 

be 


where V]_ and v 2 are now the velocity components in the direction 
of a-j_ and ct 2 , respectively. On the other hand, upon referring to 
the first of equations (l), it is still true that w n = 0, and on account 
of this it is evident that: 
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1 Sv 2 
Vl = hi BT 


( 6 ') 


If the component of velocity in the direction of the radius vector is 
expressed as a periodic function of 6, as is still permissible, then 
equation (6') immediately furnishes the means of obtaining the corre- 
sponding expression for v^. 

h. It is now easy to determine how to continue the construction of 
the flow field downstream of a given conical reference surface, Z, upon 
which the physical conditions of the flow are assumed known. Let the 
equation of the conical reference surface, Z, be given by 

cp = cp(0) (7) 


From the relationship 




Bs 

Bcp 


wherein 



one obtains 


Bs 

Bcp 



Vq ; 

vgfp - Vqj sin cp 


( 8 ) 


provided cp is expressed as a function of 
above relationship allows one to calculate 

surface E . 


9 as in equation (7). The 

— at the points of the 

Bcp 


In like manner, by use of the equation 
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and by setting 


<Hp 


+ v v 


(9) 


it is found that 




9% + cpv r = A x - cpR-L 


( 10 ) 


wherein the A^ are calculated at the points of the conical surface 
through means of the relation 


fry 


dd / z 


+ cpv T 


( 10 ') 


Finally, it is found (the intermediate steps are omitted, and just 
the final result presented) that: 


dv 9 _ A . 9* 

de ^ sin q> 


Rn 


and 


Aq = 


dv/ 


II 

_ Vm _ v 0 2 _ v 9 1 . 

^ v (p Vq) sin cp ^ 

Ap and 

A^ are the expressions: 

■ * (■ 

2 'Sq 

\r rfTl x A ^ 1 OO _ j _ 

sin cp l 

• A cos i — — — — sin 

7 R Vq &p 


v r v_ 

Vn sin cp + — 1 

V~s V 


9 


( 11 ) 


9 


( 11 ’) 
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a 3 = 




sin cp 


( 11 ") 


and their values are known on the reference conical surface, E, and thus 
so also is the value of Rg. 

By means of equation (5), therefore, one obtains 


1 


c2 sin cp c2 



sin 2 cp 



Rl = 


v v 
9 cp 


sin cp c £ 


A 1 - v cp 


>1 + 


v 0 2 ) 


cot cp 




A 2 

sin cp, 


8s 

7 r 5p 


V(pV r 

c2 


rwg 


( 12 ) 


Thus it is possible to calculate the values of R 7 at the points of 
the reference conical surface, Z. This formula is the natural extension 
of the analogous relationship already derived in reference 2 in the case 
of an irrotational flow. 

The complete solution of the problem as to how to continue constructing 
the flow field downstream of the reference surface Z is thus presented 
by the formulations given as equations (8), (ll), (12), and the additional 
equation (13), since it is also evident that 


8Vg 

8cp 


1 

sin cp 



cpR^ - Vg cos cp 


c£ 1_ 8s 

7R v e Sep 


sin cp 


(13) 
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■”5. The graph of the hodograph corresponding to an arbitrary gen- 
eral semiplane, 9 = const., is represented in figure 2. Let P be 

any point whatsoever on this hodograph, and let it signify the vec- 
tor R = R-^r - RpT. In addition let r again denote the versor that 
has the sense and direction of the radius vector in the semiplane in 
question, and let t now be the versor normal to r in the semiplane 
and oriented in the sense of an increasing cp. With these conventions 

dP 

it follows that R x — = 0 and therefore it becomes clear that the 

dp 

tangent to the hodograph at the point P is perpendicular to the above- 
defined vector R . Thus this tangent is inclined to the direction of 
r by an angle 


-1 Ro 

X = tan == ( 14) 

R 1 

If the linear element of length along the hodograph is denoted by ds 
then the absolute value of ds/dp is given by 


= ^/r-l 2 + Rp 2 

Thanks to the formulae developed in the preceding sections the 
drawing in of the tangent to the hodograph at the point P, situated on 
this hodograph, is therefore made possible, since everything is known 
about the point P. If the element of length along the hodograph is 
calculated as 


ds 

dp 


ds = l/R-j 2 + Rg 2 dcp 

then it is easy to find another neighboring point P’. When used repe- 
titious ly in the various semiplanes, this procedure allows one to determine 
the respective hodographs as the envelope of their tangent lines. 

The angle of intersection da with respect to the direction of the 
linear element of the hodograph ds is thus determined by 
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da = dcp 



R 


dRr- 


*P 1 - 


1 W~ " *2 


dR 1 

dcp 


Rq + *2 


The radius of curvature of the hodograph at the point P just 
mentioned turns out to be consequently 



When the point P' is determined in the manner that was described 
above, in consequence of having available all information about the known 
point P, it follows that Rq and Rr, are also known at the point P', 

dRq dR? 

and thus the values of and are calculable. Through means of 

dtp dcp 

equation (15) the radius of curvature R c is then determined, using 
dRq dRo 

these values of and . The direction of the principal normal 

dcp 

as defined by equation (l4) is also known, and thus the hodograph can 
be obtained in every meridional semiplane as the envelope of the respec- 
tive osculatory circles. 

6. The results obtained here for the determination of the conical 
field of flow in the case of rotational motion are employed in ah exactly 
analogous way as described in reference 2 when making a numerical appli- 
cation. In general it will be convenient to assume as the reference 
surface, Z, upon which the initial values are taken to be known, the 
conical surface of the shock wave. The required information about the 
shock wave may be obtained in first approximation by utilization of the 
hypothesis that the flow is irrotational. 


Translated by R. H. Cramer 

Cornell Aeronautical Laboratory, Inc. 

Buffalo , New York 
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Figure 1.- Coordinate and vectorial orientation in the spherical coordinate 

system. 



Figure 2.- Construction in the hodograph plane for finding velocities at 

P 1 when they are known at P. 
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